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Abstract 


We  consider  the  form  of  the  radio-frequency  (RF)  or  dielectric  input  expression  in  the  heat  equation  which 
arises  in  the  modeling  of  the  curing  of  epoxy  adhesives  in  bonding  of  composites.  We  review  two  standard 
derivations  of  a  commonly  used  expression  for  the  RF  heating  source  term.  In  this  context  we  discuss 
difficulties  involving  the  associated  inherent  polarization  assumption  and  asymptotic  behavior  of  dielectric 
parameters  as  a  function  of  the  frequency  of  the  RF  signal.  These  difficulties  cast  doubt  on  the  validity  of  the 
standard  RF  expression.  We  then  use  two  standard  polarization  models  (Debye  and  Lorentz)  to  demonstrate 
how  one  can  systematically  derive  general  RF  heating  expressions  which  do  not  suffer  the  inconsistencies 
that  arise  in  the  standard  RF  expression. 
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1  Introduction 


The  use  of  radio-frequency  (RF)  electromagnetic  energy  for  the  curing  of  adhesives  is  a  commercially 
important  process.  The  complexity  of  this  process,  one  involving  phase  transitions  and  exothermic 
chemical  reactions,  provides  a  number  of  interesting  scientific  and  technical  questions  and  challenges.  The 
application  of  interest  involves  the  adhesive  bonding  of  sheet  molding  compound  (SMC),  a  glass  fiber 
reinforced  polyester  commonly  used  in  the  automotive  industry  for  exterior  body  panels.  In  this  process,  a 
thin  layer  of  nonmagnetic  adhesive  is  placed  between  two  sheets  of  SMC.  This  adhesive  ’’sandwich’’  is 
placed  between  two  aluminum  electrodes  which  contact  the  part.  A  high  voltage  (e.g.,  1-10  KV)  electric 
field  is  applied  across  the  part  at  30  MHz.  The  field  passes  through  the  sheet  of  SMC,  which  is  relatively 
non  dielectric  (i.e.,  a  conductor),  and  into  the  dielectric  adhesive.  The  adhesive  heats  and  eventually  cures 
to  bond  the  sheets  of  SMC  together.  One  interesting  computational  challenge  is  to  model  and  predict  the 
heating  of  adhesives  so  that  precise  control  algorithms  can  be  developed  for  the  production  process.  We 
have  in  [Banks,  et.  al.,  1998]  developed  such  a  model  along  with  an  efficient  computational  simulation 
method. 

The  model  used  to  describe  the  heating  of  the  adhesive  is  a  3-dimensional  nonlinear  heat 
equation  [Bejan,  1993,  p.  16],  [Malaczynski  and  Cip inski,  1996,  p.  107] 

dT 

pcr  —  -  V  •  (k  •  VX)  i  qex  +  qrj 

where  p  is  mass  density,  cp  is  specific  heat,  pcp  is  volumetric  heat  capacity,  and  k  is  thermal  conductivity  of 
the  adhesive  which,  of  course,  may  all  be  dependent  on  T  (e.g.,  see  [Banks,  et.  al.,  1998]).  The  term  qex 
represents  the  rate  of  heating  due  to  an  internal  exothermic  reaction  and  qrj  is  the  rate  of  heat  generated 
from  conversion  of  electrical  energy  to  molecular  vibrational  energy.  (Here  the  overdots,  which  are 
standard  notation,  represent  rates,  not  actual  time  derivatives).  Our  discussions  here  focus  on  the  form  of 
the  qrj  term.  Other  aspects  of  the  model,  including  a  detailed  explanation  of  the  exothermic  term  qex, 
which  depends  on  the  degree  of  cure  and  temperature  in  a  nonlinear  manner,  are  discussed  in  the 
paper  [Banks,  et.  al.,  1998]. 

In  this  note  we  address  issues  related  to  the  RF  input  term  qrj  and  the  validity  of  a  commonly  used 
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form  of  this  term.  In  our  discussions  we  offer  possible  alternatives  to  the  usual  expression  found  in  the 
research  literature  along  with  suggestions  as  to  why  these  alternatives  might  be  more  appropriate  in  certain 
applications. 

2  Review  of  derivation  of  standard  RF  term 

In  investigating  the  term  qrf ,  one  encounters  the  following  form  which  is  widely  used  in  literature 
[Bunget  and  Popescu,  1984,  p.  17],  [Frohlich,  1958,  p.  14],  [Malaczynski,  1988,  p.  1271], 

[Malaczynski  and  Cipinski,  1996,  p.  107],  [Ross,  1982,  p.  5824],  [Scaife,  1989,  p.  13],  [White,  1973,  p.  48]: 

q  —  Lt)£ae'R  tan6EftMg.  (1) 

Here  ui  is  the  frequency  of  the  RF  bonder,  So  is  the  permittivity  of  free  space,  s'R  =  is  the  relative 
dielectric  constant  where  s  =  e'  —  is"  is  the  complex  permittivity,  tan^  =  p-  is  the  dielectric  loss  tangent, 
and  Erms  is  the  root  mean  square  strength  of  the  electric  field  being  applied  to  the  part.  Here  we  consider 
first  the  method  found  in  [Frohlich,  1958]  to  derive  (1). 

2.1  Frohlich’s  Derivation 

Much  of  the  literature  tacitly  assumes  that  dielectric  heating  is  accomplished  through  dipolar  rotation  and 

we  shall  make  this  assumption  for  the  present.  This  means  that  energy  dispersion/conversion  via 

polarization  is  of  foremost  importance.  A  first  step  in  deriving  (1)  is  understanding  the  assumptions  which 

// 

lead  to  the  definition  tan  S  =  ~j-.  Following  Frohlich,  we  consider  the  case  in  which  one  applies  a  periodic 
electric  field  E  =  Eq  cos  tot  where  Eq,  the  magnitude  of  the  electric  field,  is  independent  of  time  and  w/2tt 
is  frequency  in  cycles  per  second.  In  this  note  we  will  denote  the  electric  flux  density  (or  displacement)  by 
D  =  cq SooE  +  P  where  P  is  the  polarization  and  s ^  is  the  upper  limit  for  e  asw  — >  oo  [Frohlich,  1958,  p. 
5],  This  differs  from  much  of  the  literature  where  the  convention  is  D  =  s.qM  +  P,  but  this  change  is 
necessary  to  ensure  that  the  time  domain  polarization  formulations  are  consistent  with  the  frequency 
domain  formulations  often  found  in  literature.  Additionally,  in  some  of  the  literature,  derivations  are 
carried  out  in  the  Gaussian  system  where  So  is  not  present  in  the  definition  of  D.  Since  our  considerations 
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are  carried  out  in  the  inks  system  they  may  differ  slightly  (by  scale  factors)  from  those.  As 
in  [Frohlich,  1958],  we  assume  that  if  E  is  periodic  in  time,  P  will  also  eventually  be  periodic  in  time.  Then 
D  will  also  be  periodic  in  time  since  D  =  SoSooE  +  P .  (In  general  this  is  not  true.  For  example,  in  a  Debye 
or  Lorentz  medium,  P  is  only  asymptotic  to  E.) 

In  the  Frohlich  derivation  one  also  assumes  that  while  D  will  be  of  the  same  period  as  E,  it  will  not 
necessarily  be  in  phase  with  E.  The  rationale  is  that  matter  does  not  respond  instantaneously  to 
excitement  induced  by  electromagnetic  waves,  but  will  respond  with  some  hysteresis.  This  phenomenon 
can  be  expressed  by  the  widely  used  expression  (e.g.,  see  [Stuerga  and  Gaillard,  1996,  p.  95]) 

P(t)  =  f  K(t  —  s)E(s)ds, 

Jt—r 

which  embodies  the  assumption  that  the  polarization  at  time  t  depends  on  the  electric  field  over  an 
immediate  past  period  of  length  r.  If  we  consider  a  two  point  quadrature  (for  example,  a  trapezoid  rule  for 
the  integral  -see  [Burden  and  Faires,  1993,  p.  176])  using  the  endpoints,  from  this  polarization  assumption 
we  obtain  an  approximation  of  the  form 

P(t)  »  A'i  E(t)  +  K2E(t  -  r) 

where  K\  =  i)/i(0)  and  K2  =  |A(r).  This  means  that  P  depends  on  E  at  a  past  time.  Therefore  there  is 
a  delay  in  P  and  hence  in  D.  This  delay  is  represented  by  a  phase  shift  6.  Let  Dq  denote  the  magnitude  of 
the  electric  displacement  D.  Then 

D  = 

D  = 

where  D\  =  Do  cos  6  and  D2  =  Do  sin  6. 

For  most  dielectrics,  Do  will  be  proportional  to  Eo  and  Do/ Eo  will  depend  on  factors  such  as 
frequency,  amplitude,  temperature  and  humidity.  For  the  Frohlich  derivations,  however,  one  assumes  that 


Do  cos  (cot  —  6) 

Do  (cos  tot  cos  6  +  sin  ut  sin  6 ) 

(Do  cos  6)  cos  tot  +  (Do  sin  6)  sin  ut 

D\  cos  tot  +  D2  sin  ad  (2) 
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s'  and  s"  depend  solely  on  frequency.  We  introduce  s'  (to)  and  s"(uj)  as 


Di  —  s'  Eo 


IE  =  s"E n 


and  note  that  the  frequency  dependence  of  s'  and  s"  is  not  directly  an  issue  for  our  application  to  curing  of 
adhesives  since  we  work  at  only  one  frequency  in  our  experiments.  We  have  D\  =  s' Eq  =  Do  cos  8  and 
_t>2  -  s" Eo  =  Of)siii  8  and  consequently: 


s"  Do  sin  8 
s'  Do  cos  8 


=  tan  8. 


More  generally,  we  note  that  for  a  general  electric  field  E(t)  —  Eoelwt ,  the  formulation  from  (2)  and  (3)  is 


equivalent  to 


D  =  DeSv*-^  =  (Do  cos  8  -  Wot  sin  8)^' 


=  (s'  -  is")Eoetwt  =  sE 


where  s  =  s'  —  is"  is  the  complex  permittivity. 

We  continue  the  derivation  of  (1)  using  the  equation: 


dU  =  dQ  +  EdD 


where  dU  is  the  increase  in  total  energy  per  unit  volume,  dQ  is  the  influx  of  heat  per  unit  volume,  and 
EdD  is  the  influx  of  electromagnetic  energy  into  the  dielectric  per  unit  volume  if  the  displacement  D  is 
increased  by  dD  (see  Section  2.1.1  below  for  an  explanation  of  this  term)  possibly  due  to  a  change  in  the 
field  strength  E  at  a  constant  temperature  or  to  a  change  in  temperature  at  a  constant  E  [Frohlich,  1958, 
p.  11].  We  calculate  the  amount  of  electric  energy  transformed  into  heat  on  an  average  of  one  period.  Since 
we  assume  that  the  total  energy  is  kept  constant  and  E  is  periodic,  U  must  be  constant  and  hence  dU  =  0 


in  equation  (4);  this  yields 


dQ  =  -EdD 


and  hence  over  a  period  of  [0,  2tt/w\  we  have 


l>2Trjuj  i>2tt j u)  y*2 tt/oj  f)  T1 

dQ  =  —  EdD  =  —  E——dt. 

3  Jo  Jo  dt 
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Then  the  average  (over  one  period)  rate  of  loss  L  of  energy  from  the  electric  field  (the  heat  produced  per 
second  per  unit  volume)  is  given  by 


/»27t/u> 


r’l-K/w  /•-"/“  ol) 

L  =  ~  —  d,Q  = - /  dQ  =  —  /  E—dt. 

tJv  2tt  k  ^  Jo  dt 

We  recall  E  =  Eq  cos  Lot ,  and  from  (2)  we  have  D  =  e'Eq  cos  Lot  +  e"  Eq  sin  cot,  and  hence, 

loEq  (  Ed,  ,  ,  ,  f2*/w  „  ,  ,  \ 

— —  /  —  —  (cos2  wt)dt  +  /  e  uj  cos2  cot  dt  ) 

2^  yJo  2  at  Jo  J 


r2w/w  dD 


L  = 


r2lt/w  d  cos  tot  ..9  sin  cot, 

cos  L0t(e' — — —  +  s"——)dt 


6  2  4|2tt/uj  <"  '•(:  III  "  \ 

-2TTc”  “t.  +— £(-l 

£"E2lo 


Since  e"  =  e'tanii,  we  thus  have 


(5) 


L  =  —  s' (tan  6)  EqLO 


1  s'  i 

—So — loE q  tan  6 

2  £o 

1 


=  —cqs'j^loEq  tan  6. 


If  we  let  Erms  =  Eo/y/2  then  we  have 


q  —  L  —  £qS rLO E^m s  tan  S 


which  is  exactly  equation  (1). 


2.1.1  Origin  of  EdD  term 

We  return  to  consider  the  rationale  behind  using  the  term  EdD  as  we  did  in  section  2.1  (e.g.,  the  rationale 
behind  (4)).  We  begin  with  Maxwell’s  Equations: 


Vx£ 
V  x  H 


_dB_ 

~~dt 

dD 

~m+J 


(6) 

(7) 


V  •  D  =  p 


V  •  B  =  0 
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where  D  =  D(E)  and  B  =  B(H).  If  we  multiply  (7)  by  E,  (6)  by  H  and  subtract,  we  obtain 


E  ■  (V  x  H)  -  H  •  (V  x  E)  =  E  •  +  J)  +  H  •  (^). 

Using  the  identity  A  ■  (V  x  B)  —  B  ■  (V  x  A)  =  V  •  (B  x  A)  we  write 


.  .  ,dD  ,dB. 

-y.(E,H)  =  E.{-  +  J)  +  H.(-) 


which  can  be  rewritten  as 


d  D  dB 

[E.  ■  —  +  H  ■  — )  +  V  •  (E  XH)  +  E :  J  =  0. 


If  we  assume  that  B  =  (IqH  , 


(E-^  +  fioH-  ^+^5+^  =  0 


^ ^  m  (=) 

(a) 

where  S  =  E  x  H  is  the  Poynting  vector.  The  term  (b)  gives  us  the  rate  of  efflux  of  E-M  energy  per  unit 
volume.  The  term  (c)  gives  us  the  rate  of  conversion  of  E-M  energy  into  other  types  of  energy  per  unit 
volume.  Therefore,  by  conservation  of  energy,  (a)  must  represent  the  rate  of  change  of  energy  content  per 
unit  volume.  We  further  note  that 


rr  d-ff  _  ^0  d  2\ 


d  hqH2 
dt  2 


where  4/ro-ff2  is  the  density  of  magnetic  energy.  As  a  result  it  is  usually  suggested  that  f  E(D)dD  is  the 
density  of  electric  energy,  which  finally  suggests  that  Ed.D  is  the  change  of  energy  density  (per  unit 
volume)  connected  with  a  variation  of  D  by  dD. 


2.2  Necessity  for  alternate  derivations 

In  studying  the  derivation  above  and  the  resulting  expression  (1),  we  encountered  several  problems  with 
both  the  intuitive  behavior  of  the  RF  term  at  high  frequencies  and  with  the  implicit  assumptions  made 
about  polarization  underlying  the  derivation.  We  discuss  these  difficulties  here  to  establish  that  there  is  a 
necessity  for  alternate  derivations  and  interpretations  for  qrj. 

In  general  we  have  D  =  eoSooE  +  P  =  sE  or  in  frequency  domain  D(lo)  —  e(uj)E(u})  where  D  and  E 
denote  the  Laplace  transforms  of  D  and  E,  respectively  and  the  complex  permittivity  is 


e(w)  =  £re(u>)  +  i£im(u>)  —  £'  —  i£n .  We  note  that  (1)  can  be  written  as  qrf  =  ujP'  E^MS.  As  w  — -  0  we 
approach  a  static  field  for  which  there  is  no  dielectric  loss  and  the  permittivity  approaches  a  constant 
denoted  £q£s  where  £s  is  the  static  dielectric  constant.  As  a  result,  we  expect  to  see 

e"(ui)  — *  0  as  w  — >■  0 


(8) 


e'(u j)  — »■  £o£s  as  u  — >  0. 

Continuing  our  heuristic  intuition-based  arguments,  we  consider  e'  and  e"  as  w  — ^  <30.  Clearly  we  do  not 
want  e"  to  become  unbounded  as  ui  becomes  unbounded  and  so  we  require  £n  to  approach  some  constant. 
However,  if  e"  approaches  a  nonzero  constant  then  qrf  will  become  unbounded  as  w  becomes  unbounded, 
which  is  undesirable  behavior  since  we  expect  a  finite  bound  on  qrj  for  a  given  material.  Hence,  the 
limiting  constant  for  e"  must  be  zero  or  e"  — »  0  as  to  — >■  oo.  All  materials  have  inertia  and  hence  as  w  ^  oo, 
the  dielectric  material  does  not  move  and  the  RF  electromagnetic  signal  must  move  at  the  speed  of  light. 
Therefore,  by  definition  of  £oo  we  expect 

£"(to)  — >  0  as  to  — >■  oo 

(9) 


e'(w)  — >  £o£oo  as  ui  — >■  oo. 

Thus,  if  (1)  is  a  correct  expression  we  intuitively  expect  the  limiting  behavior  for  £',£"  given  in  (8),  (9).  As 
we  shall  see  below,  this  can  lead  to  inconsistencies. 

A  slightly  different  approach  to  required  asymptotic  behavior  of  dielectric  parameters  involves 
consideration  of  general  polarization  models  in  the  frequency  domain  of  the  form 


P(S)=£0^T(S), 

ns) 

where  P  and  E  are  Laplace  transforms  of  P  and  E,  respectively,  and  p(s),q(s)  are  polynomials  in  s. 
Letting  s  =  iu,  from  D  =  £o£oo E  +  Fwe  have 


D  M  = 


£q£c 


£0p(iu 


SqSco  + 


E(cv) 
£o  q{iu)p{iu) 


q(iu 


E(lo), 


kNI2 

where  z  and  \z\  denote  the  usual  complex  conjugate  and  amplitude,  respectively,  of  z.  If 

(degree  of  p)<(degree  of  q)  then  D  — *■  £o£c yoE  as  to  — >■  oo  and  we  have  what  is  sometimes  defined  as  a  stable 

material  polarization  model.  We  note  that  a  model  must  be  stable  in  order  for  it  to  satisfy  (8)  and  (9). 
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2.2.1  Polarization  assumptions 


In  this  section  we  describe  polarization  assumptions  implicit  in  the  derivation  of  Frohlich  to  obtain  (1). 
There  are  two  ways  to  view  the  polarization.  The  first  is  to  write  P  as  a  function  of  E  (this  will  result  in 
(11)  below),  while  the  second  is  to  write  it  as  a  function  of  both  E  and  E  (this  results  in  (12)  below).  We 
begin  again  with  D  —  SoSooE  +  P  and  recall  that  we  are  working  with  the  periodic  case  in  which 
E  —  Eq  cos uit  and  D  =  Do  cos {wt  —  8).  We  have 


Dq  COS  (wt  —  8 )  =  SqSooEq  cos  uit  +  P 


or 


P  =  Dq  cos  (wt  —  8)  —  soSooEq  cos  uit. 


We  assume  Do  is  proportional  to  Eo  (i.e„.  Do  --  0£o  where  ^  0  may  depend  on  ui)  and  substitute  into 
the  above  equation.  In  the  case  where  we  write  P  as  a  function  of  E  we  have 


P  =  0Eo  cos{uit  —  8)  —  SoSooEo  cos  uit 

=  0Eo(cos  uit  cos  8  +  sin  uit  sin  8 )  —  SoSooEo  cos  ut  (10) 

=  0E  cos  8  +  P^Eq  —  Eq  cos2  Lot  sin  <5  —  eoSooE 
=  pB  cos  6  +  f3yj Eq  —  E2  sin  S  —  Soe^E 
=  E(p  cos  6  -  eo£oo)  +  P\j Eq  -  E2  sin  6. 

Thus  we  have 

P  =  E(P  cos  6  —  eoSoo)  +  ft\j Eq  —  E2  sin  6,  (11) 

from  which  we  immediately  observe  that  P(—E)  =£  — P(E )  unless  8  =  0  or  0  =  0.  If  6  =  0  (i.e.,  D  =  s /. ) . 
we  have  P  =  E{0  —  So^oo)  —  x£  which  is  the  standard  polarization  law  for  a  linear  (sometimes  also  called 
ideal)  medium  [Reitz,  et.al.,  1993,  p.  110]  and  no  apparent  paradoxes  arise.  However,  if  we  reverse  the  field 
E,  in  a  dipolar  medium  we  must  reverse  the  polarization  P ,  which  implies  P(—E)  =  —P(E)  and  P  =  P{E) 
must  be  odd  in  E.  As  we  have  already  observed  from  (11),  for  8  ^  0,  0  ^  0,  this  is  not  true  of  the 
polarization  implicit  in  the  Frohlich  derivation  of  (1).  Therefore  (11)  is  not  a  reasonable  or  valid 
polarization  constitutive  relation  for  any  dipolar  material. 
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On  the  other  hand,  if  we  begin  with  (10)  and  proceed  in  a  slightly  different  manner,  we  also  obtain 


inconsistencies.  We  have 


=  /3Eo( cos  Lot  cos  8  +  sin  tot  sin  8)  —  sqSooEq  cos  Lot 

=  E(/3  cos  8  —  sqSoo)  —  - — -(—  cos  Lot)  sin  8 


E{j3  cos  8  —  SoSoo )  —  —  ( —  /.'c  cos  Lot)  sin  8 
to  dt 

E(Q  cos  8  —  SoSoo )  —  - — -r~  sin  8. 

lo  at 


Thus  this  argument  leads  to 


P  =  aE  +  bE 


where  a  =  f3  cos  8  —  So^oo  and  b  =  —  £■  sin  8.  The  second  potential  inconsistency  can  be  illustrated  by 
studying  the  form  of  s  implicit  in  (12).  We  take  the  Laplace  transform  of  (12)  to  obtain 


P  =  (a  +  bs)E 


where  a  =  /3  cos  8  —  SqEoo  and  b  =  —  sin  8.  Hence 


D{lo)  —  (so£oo  +  a  +  ibuj)E 

=  [(eo£  CO  H-  a)  +  i{bw)]  E 

0uj  1  - 
=  (f3  cos  8 +  i(—  ■ — sin  8)  E 

LO 


—  f3  cos  S  —i  (3  sin  S  E 


which  relates  /?  to  the  expressions  for  s'  and  e"  defined  in  Section  2.1.  Recalling  the  expected  asymptotic 
behavior  of  s'  and  e"  in  Section  2.2  leads  to  inconsistencies  for  f3  in  light  of  (13).  To  see  this,  we  note  that 
in  the  Frohlich  derivation  it  is  not  explicitly  stated  whether  or  not  the  phase  shift  8  is  frequency  dependent. 
From  (3),  it  is  unlikely  but  possible  that  tan  8  is  independent  of  lo  even  though  s'  and  e"  generally  are  not. 
If  8  is  not  frequency  dependent  (i.e.,  8  ^  8(lo))  we  see  from  (13)  that  in  order  for  (8)  to  hold  we  must  have 
simultaneously  (3  — >  and  /?  — >■  0,  as  lo  — >■  0,  which  leads  to  a  contradiction.  A  similar  contradiction 

arises  for  (9).  Thus  if  8  ^  8(w ),  (8)  and  (9)  are  violated.  In  fact,  8  should  depend  on  frequency  for  a 


number  of  reasons,  among  them  that  the  delay  time  should  increase  as  frequency  increases.  If  we  assume 
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6  =  S(lj)  then  the  asymptotic  behavior  must  be  of  a  certain  form  in  order  to  prevent  apparent 
contradictions.  This  will  be  detailed  below. 

In  a  number  of  cases,  permittivities  based  on  other  polarization  models  do  meet  the  requirements  of 
(8)  and  (9).  We  illustrate  this  with  examples  involving  two  common  polarization  models.  We  first  consider 
the  Debye  model  [Reitz,  et.al.,  1993,  p.  514]: 


p  _|_  Ip  _  go(gs  gc°)p 

T  T 

where  r  is  the  relaxation  time,  es  is  the  static  dielectric  constant,  and  s^  is  the  limiting  value  of  the 
permittivity  as  the  frequency  becomes  unbounded.  We  take  the  Laplace  transform  to  obtain 


(«  +  -  )p  =  £o(£s  £oo)it 

T  T 


Hence 


d=[£0£oo+£o^  e™y]E. 

1  +  TS 


Taking  s  —  ico  we  have 


D(cv)  = 


g0  goo  + 


g0gc 


go(gs  goo) 


E 


1  +  iuir 
go(gs  -  goo)(l  -  Wt) 


1  +  lo2t2 


E 


go(gs  goo)>\  .  ( g0 ^T~(6s  €qo ) 

£o£°°  +  1  +lj2t2'  )  1  (  l+w2r2 

-  -  ^ 


E. 


(14) 


As  w  — >■  0,  it  is  clear  that  e'  — »■  £q ss  and  e"  — ?  0.  As  w  becomes  unbounded  we  see  that  s'  — *  £o£oo  and 
s"  0,  which  is  exactly  the  behavior  required  by  (8)  and  (9). 

We  next  consider  the  Lorentz  model  which  is  described  in  detail  in  [Reitz,  et.al.,  1993,  pp.  495-499] 
given  by 

P  4 — P  T  ojf)  P  =  StjujpE. 

T 


Here  r  is  again  the  relaxation  time,  loo  is  the  natural  frequency  of  the  undamped  oscillator,  and 

Up  =  Ne/som  where  N  is  the  number  of  charges  per  unit  volume,  e  is  the  displaced  charge  and  m.  is  mass 

of  the  particle  in  the  damped  harmonic  oscillator  paradigm  that  is  the  basis  of  the  Lorentz  model.  Again, 
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we  take  the  Laplace  transform  and  obtain 


“ P  H — sP  +  LO qP  —  £qLO“E 


and  taking  s  —  ito  we  have 


P  = 


£0  TLO2 

r(u>Q  —  LO2)  +  iu> 


E. 


Using  again  D  =  £g£t x>E  +  P  we  find 


D 


£o£<x>  + 


£q£c o  + 


£otu)£ 


(Wq  -  LO2)  + 


£c(r(wo  -  w2)  -  iw)rw2 
r2(u>2  -w2)2  +w2 


M 


(  £0T^-uP)^\  ,( 

^£0£oo  +  r2(w2  _  w2)2  +w2  J  *  ^r2( 

e' 

As  w  ^  oo  we  see  that  s'  — *  £o£oo  and  s"  — *  0.  Additionally,  as  w  — *  0,  e"  — >  0,  but 

c  t-2,  ,2,  ,2 

/  _  ,  £0T  WqOIj, 

£  —  £q£oo  I - ^ — 4 — 

r"u>0 

and  it  is  not  clear  that  this  is  equal  to  £q£s.  However,  if  we  let  lo2  =  Wq(£s  —  £rc)  [Scaife,  1989,  p.  59],  then 
s'  — >  £o£s  as  required  in  (8)  and  the  Lorentz  model  satisfies  (8)  and  (9). 

We  turn  our  attention  to  the  form  of  qr  f  in  ( 1 ) .  Clearly  qr  f  — *■  0  as  lo  — >  0  since  we  know  s"  — >■  0  as 
lo  — >  0.  As  w  becomes  unbounded,  however,  the  behavior  of  qrj  is  not  as  clear.  Since  we  want  qrj  to 
approach  some  constant  as  u  becomes  unbounded,  we  want  s"  — *■  0  faster  than  lo  — >  oo.  In  order  for  this  to 
hold,  s"  must  behave  like  i+a  for  some  a  >  0.  Similarly,  tan  A  must  behave  like  i+a  ■  If  this  is  true,  as  lo 
becomes  unbounded  qrj  will  behave  like  ^  which  approaches  0  for  a  >  0  and  approaches  a  constant  for 
a  =  0.  If  this  is  not  true,  qrj  will  become  unbounded  as  lo  becomes  unbounded,  and  we  have  an  objection 


SgLOTLO 


V 


LO2  -  LO2)2  +  LO2 


E. 


(15) 


to  (1). 


We  note  that  stability  is  not  enough  to  guarantee  that  qrj  will  approach  some  constant  as  lo  becomes 
unbounded.  Stability  guarantees  only  that  s"  approaches  0,  giving  no  information  on  the  rate  at  which  s" 
approaches  0  relative  to  the  rate  at  which  lo  becomes  unbounded.  Therefore,  we  have  a  stronger  restriction 
on  the  materials  whose  behavior  can  be  approximated  by  (1).  In  Section  3  we  derive  alternatives  to  (1) 
which  satisfy  this  added  restriction. 
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2.3  Alternate  derivation  of  standard  RF  term 


Since  we  encounter  difficulties  with  the  implicit  polarization  assumption  in  the  derivation  in  Section  2.1,  we 
present  a  derivation  of  (1)  with  arguments  that  do  not  involve  the  polarization.  This  derivation,  based  on 
the  report  [Lord,  1983]  and  arguments  in  [White,  1973],  relies  on  Joule’s  heating  law  along  with  a  number 
of  other  somewhat  questionable  tacit  assumptions.  Again  one  obtains  the  standard  RF  heating  rate  term 
(1),  but  from  energy  considerations.  Joule’s  heating  law  [Elliot,  1993,  p.  489],  under  assumptions  of  no 
conduction,  radiation,  or  other  losses,  i.e.,  100%  efficiency,  is  given  by 

f  =  ^  (16) 

where  U  is  internal  energy,  a  represents  conductivity,  and  E  represents  the  instantaneous  electric  field  in 
the  dielectric.  We  note  that  the  dielectric  must  satisfy  Ohm’s  Law  (J  =  crE ,  a  standard  assumption  for 
linear  isotropic  materials  [Reitz,  et.al.,  1993,  p.  167])  in  order  for  this  form  of  Joule’s  law  to  be  valid.  From 
(16)  we  have  over  any  time  period  [fg,  t\ 


r  du  _  r 

Jto  dt  Jt0 


=  /  crE~di 


or,  assuming  a  constant, 


AU  =  a  E2dt 


where  AU  -  U(t )  —  U(to). 


If  we  further  assume  a  periodic  field  E  =  Eq  cos  uit,  then  we  have 


AU  =  uEq  I  cos"  us  ds 


=  a Eq  f  — [1  +  cos  2uis]d,s 
Jtn  2 


narl  ,  1 


=  o-Eo[-s  +  —  sin2ws]<o 


=  aEl  (77O  -  to)  +  -7—  sin  2 cut  -  —  sin  2 wt0) 


where  the  second  and  third  terms  in  the  last  expression  are  very  small  (and  hence  neglected)  because  RF 
heating  involves  w  of  approximately  106  —  10'  Hz.  If  we  assume  to  =  0,  we  have 


AU  =  -aEh. 
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We  next  consider  the  conductivity  a  and  argue  that  it  is  given  by  cr  =  usos'R  tan  6.  We  note  that  if  (1)  is  to 
hold,  this  must  be  the  case. 

The  derivation  is  suggested  by  White  [1973]  who  notes  that  in  the  literature  conductivity  is 
commonly  written  a  =  us" .  Using  this  along  with  the  definition  of  s"  given  in  Section  2.1  we  see  that 


a  =  us1  tan  8 

s1 

=  us  o — tan  8 
£o 

=  usos'fttanS. 


With  this  form  for  a  and  the  expression  for  we  see  that  the  rate  of  absorption  L  can  be  written 

AU  1  2  1  .  ., 

L  —  —j—  —  -aUg  —  -u5asRE^tan8. 


Setting  Erms  —  Ssl/v/S,  we  obtain 


q  —  L  —  usqs RERM s  tan  8. 


We  see  that  this  derivation  yields  exactly  the  same  expression  as  the  Frohlich  derivation.  Although  this 
derivation  does  not  depend  explicitly  on  the  polarization  assumptions  found  in  our  earlier  derivation,  they 
are  inherent  in  our  definitions  of  a  and  therefore  this  derivation  leads  to  the  same  conceptual  difficulties  as 
the  previous  one.  Moreover,  the  assumptions  (implicit  as  well  as  explicit)  inherent  in  the  arguments  (e.g., 
constant  conductivity,  ohmic  rather  than  dielectric  heating,  etc.)  really  are  not  viable  for  most  dielectrics. 


3  Alternate  RF  expressions 

At  this  point  we  briefly  summarize  the  difference  between  ionic  and  dipolar  conduction  as  presented  in 
[White,  1973].  In  general  adhesives,  ions  naturally  exist.  Concentration  is  dependent  on  chemistry, 
presence  of  water  and  other  impurities  or  contamination.  These  ions  are  electrically  charged  and  as  a  result 
are  attracted  or  repulsed  by  electric  fields.  The  ions  then  collide  many  times  (depending  on  the  number  of 
random  ions,  density  of  material  and  the  frequency  of  the  field)  with  other  nn  ionized  molecules  which 
provides  a  two  step  conversion  process  in  which  electric  field  energy  is  converted  into  ordered  kinetic  energy 
and  then  into  disordered  kinetic  energy  where  it  is  appropriate  to  regard  it  as  heat.  Dipolar  conduction,  on 
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the  other  hand  takes  place  due  to  the  dipolar  rotation  of  the  molecules  in  the  material.  Here  there  is  a 
conversion  of  energy  from  electric  field  energy  to  potential  energy  as  the  electric  field  builds  up  and  then  to 
stored  random  kinetic  or  thermal  energy  in  the  material  as  the  field  dies  down  [White,  1973,  p.  47].  We 
feel  that  the  epoxy  system  which  we  are  modeling  has  very  low  ionic  content.  Hence,  for  our  derivations 
below  we  will  assume  all  heating  is  dipolar,  but  in  reality,  depending  on  the  material,  there  may  be  other 
ionic  contributions. 

Since  we  have  shown  that  different  derivations  of  (1)  depend  implicitly  on  a  polarization  model  with 
some  inconsistencies,  we  now  derive  alternate  forms  of  the  qr  f  expression  based  on  the  Debye  and  Lorentz 
polarization  models  often  assumed  in  the  electromagnetic  materials  literature.  We  derive  first  a 
Debye-based  heating  expression  and  discuss  its  asymptotic  behavior  and  then  do  the  same  with  a 
Lorentz-based  expression.  We  do  this  to  illustrate  the  considerations  that  must  be  incorporated  into  any 
serious  effort  to  derive  general  RF  expressions. 


3.1  Derivation  of  Debye-based  RF  expression 


We  wish  to  derive  an  RF  heating  term  based  on  a  Debye  polarization  model.  Recall  (see  Section  2.2.1)  that 
the  Debye  model  is  represented  by  the  equation 


1 

P  +  -P  = 


If  we  let  a  =  -  and  b  =  £°°)  we  have 

T  T  5 


(17) 


P  +  a.P  =  bE. 


(18) 


Again  we  assume  E  =  EocosLot  where  lo/2it  is  the  frequency  in  cycles  per  second.  Substituting  into  (18) 
we  obtain 

P  +  aP  =  bEo  cos  Lot 


which  yields  the  solution 

P(t)  =  Poe~at  +  bEoe~at  (  eas  cos  los  d.s. 

Jo 

Substituting  (19)  into  the  equation  D  =  SqSooE  +  P,  we  find 


D  =  £o£oo  Eo  cos  cot  +  Poe  ai  +  bE$e  at  f  eas  cos  los  d.s. 

Jo 


(19) 


(20) 
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As  in  the  previous  derivations  (see  (5))  we  have 


u  f2*r*  „ dD  , 

—  /  E——di. 

2  tt  J0  dt 


So  from  (20)  we  obtain 


u  f-w/w 
—  Eo  cos 

■E*  Jo 

/2^£o e 

Er  °  Jo  2 

^E2  [ 

2tt  0  Eo  Jo 


O'  ft 

I  EocosLot—  eoSooEo  cos u>t  +  Poe~at  +  bEoe~ai  /  eas  cos usds  dt 

c  O'  L  Jo  . 

»i  '  !!rLl(“s2“‘)  +  ^cos“4f‘",)+l,cos“‘l[e"'i  *“■' 


eas  cos  us  ds  dt 


2tt  j  U)  p2'Kfw  ^  rt 

e~at  cos  ut  dt  +  b  /  cos  cut—  /  eas  cos  us  ds  dt 


Carrying  out  the  remaining  integration,  we  find 


L  =  —  En 


2P0  (e-2lra/“  -  1)  b(uir(a2  +  u2)  -  a3e-2™/“  +  a3) 


2? r  u  \  Eo  a2 +u2 


(a2  +  w2)2 


and  if  we  let  Erms  —  ^=,  we  have  an  alternate  RF  term: 


q  —  L  —  u 


2P0  (e~2l,alw  -  1)  b(uir(a2  +  u2)  -  a3e~2lra^  +  a3) 


ttEo  a 2  +  u2 


tt  (a2  +  u2)2 


—  i  Emrl  h  —  g°(.£'-£") 


where  a  =  -  and  b  = 


We  study  the  behavior  of  this  term  which  can  be  rewritten 


q  =  uTE£ 


where  T  represents  the  bracketed  term  in  equation  (21).  This  has  the  same  form  as  (1),  but  the  constants 
must  be  interpreted  differently  (in  (1)  T  is  given  by  eo£^tan<5).  It  is  easily  seen  that  as  u  approaches  0, 
r  — >■  +  £o(£s  —  which  is  constant  and  hence  q  — »  0,  which  is  the  same  behavior  found  for  (1). 

As  u  approaches  infinity,  however,  we  see  much  different  behavior.  We  expand  T  in  (21)  to  obtain 


a2P0  e-2™/“  1  bua2 

irEo  a2  +  u2  a2  +  u2  +  (a2  +  u2)2 

bu3  ba3e~2waE  |  a3 

+  (a2  +  u2)2  ~  7 r(a2  +w2)2  +  tt(o2  +  u2)2  ' 

v. - v - " - - v - - >  " - - v - 

id)  C)  (/) 

Clearly  terms  (a),  (b),  (e),  and  (f)  go  to  0  as  u  — >■  oo,  and  by  applying  L’Hopital’s  rule  to  the  remaining 


terms  we  see  that  they  too  approach  0  as  u  approaches  infinity.  If  we  relate  T  to  e"  in  the  standard  term 
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and  recall  our  discussion  of  the  desired  behavior  of  e"  in  Section  2.2.1  where  we  let  £im  =  e" .  this  is  exactly 
what  we  expect  to  see.  If  we  now  look  at  the  entire  term  given  in  (21),  using  the  same  methods  for  taking 
the  limits  as  above,  we  find  that  as  oj  — >  oo, 

,  £0  (&S  &co) 

q  — -  6  =  - -. 

T 

That  is,  as  the  frequency  becomes  very  large,  the  rate  of  heating  will  approach  a  finite  constant  value,  as 
discussed  in  Section  2.2.1. 


3.2  Derivation  of  Lorentz-based  RF  expression 

Another  common  polarization  model  is  the  Lorentz  model.  Following  the  same  steps  used  Section  3.1,  we 
derive  an  RF  term  based  on  this  model.  Recall  that  the  Lorentz  model  is  represented  by  the  equation 

P  +  -P  +lo2P  =£0  uj2E.  (22) 

T  y 

If  we  let  a  =  j,b  =  lo o,  and  c  =  £o co2,  we  have 

P  +  aP  +  b2P  =  cE.  (23) 

We  assume  E  —  Eo  cos  cot  and  substitute  into  (23)  to  obtain 

P  +  aP  +  b2P  —  cE0  cos  u>t. 


After  applying  the  variation  of  constants  formula,  we  find  a  solution  P  to  (23) 

P0(eAl<A2  -  eA2lAi)  Pi(eA2*  -  eAl*) 


P(t)  = 


7 


7 


cE0 

7 


ui  sin  ad  +  A2eA2t  —  A2  cos  uit  Ai  cos  it it  —  to  sin  cot  —  Aie 


A 1 1 


\2+lo< 


A  2+w2 


(24) 


where  7  =  \/ a2  —  4 b2 ,  Ai  =  ~a~1 ,  and  A2  =  ~a^1 .  We  note  that  Ai  and  A2  are  the  eigenvalues  of  the 
system  of  first  order  equations  equivalent  to  (23)  (which  is  a  damped  harmonic  oscillator)  and  since  this 
represents  an  asymptotically  stable  dynamical  system  [Brauer  and  Nohel,  1969,  p.  151],  we  have 
Re(A,;)<  0,  for  i  =  1,2.  Following  the  previous  arguments,  we  wish  to  find  L  given  by 


L  = 


2? r 


f2,r/“  dD 

E~E7dt 
t  at 
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where  D  —  SoSooE  +  P  and  P  is  given  by  (24).  After  some  simple  integration,  this  reduces  to 


loE0  OP 

L  =  — —  /  cos  wt-^-dt. 

Jo 


2  to  Jo  dt 

Carrying  out  the  remaining  integration  and  letting  Erms  —  Jip/\/2,  after  some  straightforward  but  tedious 
calculations,  we  obtain  a  form  for  the  RF  term: 

1  ((\iMPo  A  P  cA2  \  A:(c2*A’^  1) 


q  —  L  =  ui 


ttj  \\  Eo  Eo  (Af  +lo2))  (Aj+w2) 

A •  A2 Po  A2Px  cA2  \  A2(e2'A»/"  -  1) 


Eo  Eo  {A2+lu2)J  (A2+w2) 

1  1 


At  +  w2  A?  +  < 


7712 

^ RMS • 


(25) 


We  investigate  the  behavior  of  this  expression,  which  can  again  be  written  in  the  form 


q  =  uiTEj 


RMS 


where  T  is  the  bracketed  term  in  (25).  We  first  consider  the  behavior  as  to  — >■  0.  We  note  that  as  u 
approaches  0,  e2,rAl/,“  — *  0  since,  as  stated  previously,  Re(Ai)  is  less  than  0.  The  same  holds  for  the 
exponential  involving  A2.  From  this,  it  is  clear  that  as  u  — >■  0, 


1 

r  — 

7T7 

which  is  constant  and  again  we  see  that  q  — >■  0  as  ui  — >■  0.  Applying  the  same  arguments  as  in  Section  3.1,  it 
is  easy  to  see  that  q  — >■  0  as  to  — >■  oo  which  is  the  same  behavior  found  in  the  Debye-based  term.  To  see  this 
note  that  the  terms 


AiA2Po  —  AiPi  \  1  f  — Ai  A2Po  +  A2Pi  A  1 
Wo  +c)  W  +  {  Wo  c)  w 


CTTLO2  C7TCJ2 

Af+W  and  A2  +  U)2 

each  approach  c?r  as  w  — ^  oo  while  all  of  the  other  terms  approach  0,  so  we  find  q  — »■  (c7r  —  cit)Ef>MS  =  0  as 
lo  — ^  oo.  Again,  this  limiting  behavior  is  more  intuitively  satisfying  than  the  behavior  of  the  standard  RF 
term  given  by  (1). 

Most  dielectric  materials  are  quite  complicated  in  their  polarization  mechanisms  and  the  related 
constitutive  laws  usually  cannot  be  expressed  simply  in  terms  of  a  simple  Debye  or  Lorentz  law.  However, 
it  is  quite  common  to  model  polarization  using  combinations  of  multiples  of  these  (and  other)  simple 
polarization  relationships  [Albanese,  1997],  [Oughstun,  et.al.,  1993].  Adhesives  are  no  less  complex  with 
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respect  to  polarization.  Almost  all  adhesives  (two  component  epoxies,  water  based  urethanes  as  well  as 
more  general  two  component  urethanes,  acrylics,  etc.)  involve  liquid  as  well  as  solid  components.  Moreover, 
during  the  adhesive  curing  process,  transitions  from  liquid  to  solid  phases  occur,  further  complicating  the 
derivation  of  any  polarization  based  RF  heating  law.  Thus  while  the  above  examples  illustrate  how  one 
might  go  about  deriving  an  RF  term  from  a  given  polarization  law  and  checking  it  for  consistency,  it  is 
most  unlikely  that  either  of  the  specific  examples  (Debye  or  Lorentz),  used  above  for  illustrative  purposes, 
will  provide  an  adequate  basis  for  deriving  RF  terms  for  many  adhesives  used  in  practice.  Our  current 
research  efforts  (both  theoretical  and  experimental)  are  focused  precisely  on  this  topic. 

4  Summary  and  implications 

We  have  derived  several  expressions  for  the  term  qrj .  We  see  that  all  of  these  expressions  are  of  the  form 
qrj  =  wTE2rms  where  the  definition  of  T  depends  on  the  polarization  model  used  in  the  derivation.  More 
importantly,  we  have  shown  that  the  standard  term  (1)  widely  found  in  the  electromagnetic  heating 
literature  corresponds  to  a  polarization  model  with  inherent  inconsistencies.  We  have  thus  demonstrated 
the  need  for  alternate  heating  expressions.  We  follow  the  same  steps  used  in  the  usual  derivation  to  derive 
alternate  expressions  based  on  two  well  known  polarization  models  -  Debye  and  Lorentz.  Based  on  our 
analysis,  we  suggest  that  while  the  Debye  and  Lorentz  models  are  too  simple  to  approximate  the 
polarization  behavior  of  most  adhesives,  they  at  least  result  in  RF  heating  terms  that  possess  desired 
asymptotic  behavior  properties  lacking  in  the  standard  RF  expression  found  in  the  literature.  Moreover, 
these  examples  illustrate  an  approach  that  can  be  followed  to  derive  RF  heating  expressions  based  on  a 
given  polarization  model. 

Arguments  for  the  use  of  (1)  due  to  the  simple  fact  that  s'  and  e"  can  be  found  easily  through  experi¬ 
mentation  are  not  persuasive  since  the  expression  is  not  based  on  physically  realistic  material  assumptions. 
We  can  relate  s'  and  s"  to  the  Debye  model  using  the  form  of  (14)  and  to  the  Lorentz  model  using  (15).  For 
example,  if  we  know  something  about  r  in  the  Debye  model,  we  can  use  these  relationships  to  determine  ss 
and  £ oo  and  hence  obtain  our  equation  (21)  in  terms  of  s'  and  tan  6.  Similarly,  if  we  have  information  about 
two  of  the  parameters  in  the  Lorentz  medium,  we  can  write  (25)  in  a  similar  form.  Therefore,  we  see  that 
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not  only  do  these  terms  model  the  desired  asymptotic  behavior  of  RF  heating  more  accurately,  but  with 
some  information  about  one  or  two  of  the  terms,  they  can  be  estimated  using  the  same  measured  quantities 
used  in  the  standard  RF  term  (1).  There  is,  of  course,  a  remaining  challenge  in  deriving  RF  heating  terms 
based  on  more  realistic  polarization  assumptions  for  adhesives  while  at  the  same  time  obtaining  RF  terms 
parameterized  in  a  manner  consistent  with  the  reliable  estimation  of  parameters  from  experimental  data. 
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